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Abstract 


Non-covalent functionalization via physisorption of organic molecules provides a scalable ap¬ 
proach for modifying the electronic structure of graphene while preserving its excellent carrier 
mobilities. Here we investigated the physisorption of long-chain acenes, namely, hexacene and 
its fluorinated derivative perfluorohexacene, on bilayer graphene for tunable graphene devices 
using first principles methods. We find fhaf fhe adsorption of fhese molecules leads fo fhe for- 
mafion of localized sfafes in fhe elecfronic sfrucfure of graphene close fo ifs Fermi level, which 
could be readily funed by an exfernal elecfric field in fhe range of ± 3 eV/nm. The elecfric field 
nol only creafes a variable band gap as large as 250 meV in bilayer graphene, buf also sfrongly 
influences fhe charge redisfribufion wifhin fhe molecule-graphene sysfem. This charge redis- 
fribufion is found fo be weak enough nol fo induce slrong surface doping, buf slrong enough fo 
help preserve fhe elecfronic slates near fhe Dirac poinf of graphene. Our resulls furlher high- 
lighf graphene’s polenlial for selecfive chemical sensing of physisorbed molecules under fhe 
exfernal elecfric fields. 
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Introduction 


Graphene—a planar layer of earbon atoms arranged in a hexagonal lattiee—exhibits a linear elee- 
tronie dispersion with the valenee and eonduetion bands touehing at the Dirae point. ^ As a result 
of this unique eleetronie strueture, pristine graphene demonstrates an ultrahigh eharge earrier mo¬ 
bility in exeess of 200,000 em^V^s'^ whieh ean be exploited for novel, highly energy effieient 
eleetronie deviees.^“^ However, the development of graphene-based eleetronie deviees is primar¬ 
ily hindered by the absenee of an intrinsie band gap in its eleetronie strueture.^ Although various 
approaehes for tailoring the eleetronie strueture of graphene have been pursued in reeent years, 
ereating a signifieant band gap while maintaining large eharge earrier mobilities in graphene re¬ 
mains a formidable ehallenge. 

One strategy to modify the eleetronie strueture of graphene is to utlilize quantum eonfinement 
effeets inherent in low dimensional struetures sueh as quasi one-dimensional graphene nanorib¬ 
bons. While this strategy ean effeetively induee the band gap, it also suffers from earrier seat- 
tering due to edge imperfeetions. Another route is the ehemieal funetionalization of graphene 
where the addition of eovalent bonds to graphene (for example, via hydrogenation and fluorination) 
ehanges the hybridization of earbon atoms from sp^ io sp ^While sueh eovalent funetionaliza¬ 
tion suoeessfully alters the eleetronie properties of graphene, it also leads to a severe degradation 
of its transport properties. Reeent demonstrations of heterostruetures of graphene with other 2D 
materials (for example, boron nitride and transition metal diehaleogenides) also provide a possible 
option, but a eonsistent produetion of graphene heterostrueture deviees is diffieult to eontrol on 
large seale. 

Among the various approaehes being pursued to modify the eleetronie strueture of graphene, 
non-eovalent funetionalization via physisorption of organie moleeules offers an interesting path¬ 
way. This approaeh relies on eonserving the integrity of the 5p^-bonded earbon lattiee and thus 
preserves the linear dispersion of eleetrons near the Dirae point. Moreover, the produetion of 
deviees made of graphene with physisorbed moleeules ean be readily assisted by moleeular self- 
assembly and ean therefore be expeeted to be sealable. Reeent studies have suggested that in 
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graphene physisorbed with small molecules such as NO2 and NH3, the application of a transverse 
external electric field can further enhance the tunability of the electronic structure of graphene 
by affecting the charge redistribution. Among the organic compounds that are amenable to 
physisorption on graphene, aromatic molecules are of particular interest. The face-centered 
parallel stacking of aromatic molecules on graphene surface can lead to a stable hybrid system 
via van der Waals (vdW) interactions, while the enhanced n-n electron interaction is expected 
to influence the electronic structure of graphene.^ Moreover, addition of functional groups with 
high electron or hole affinity to the aromatic molecules has been suggested as an effective ap¬ 
proach to induce strong charge doping in graphene. 22,27-30 allow a vertical integration 

of graphene with physisorbed organic molecules with tunable transport characteristics such as 
charge injection barriers. However, recent reports indicate that a strong surface charge doping of 
graphene by molecules or electric field can cause a significant shift of the Fermi level into the va¬ 
lence or conduction band, and often lead to a severe deformation of the n bands of graphene. 

As a result, the charge carrier mobility of graphene degrades, thus limiting the switching capa¬ 
bility of graphene-based semiconducting devices. Therefore, identification of suitable organic 
molecules for non-covalent functionalization of graphene still remains an open challenge for con¬ 
trollable modification of its electronic properties. 

Hexacene belongs to the group of acenes, the aromatic compounds formed by linear fusion 
of benzene rings (C4n+2H2n+4)- Long-chain acenes possess low-lying molecular orbitals that are 
expected to hybridize with n electrons of graphene and thus influence its electronic structure. 
Furthermore, it has been shown that the edges of hexacene can be readily functionalized with 
chemical groups with widely varying electron or hole affinity. Recently, Watanabe et al. reported 
a successful way to synthesize hexacene that can remain stable up to 300 °C in dark conditions. 
Moreover, organic field effect transistors (OFET) devices made of hexacene have demonstrated 
a highest charge carrier mobility ever reported for organic semiconductors.^^ These observations 
suggest that hexacene can be an attractive candidate for a stable physisorption on graphene. A 
good fundmental understanding of the electronic interactions between hexacene and graphene is 
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Figure 1: A schematic illustrating (a) single- and (b) dual-molecular adsorption of hexacene and perfluoro¬ 
hexacene on bilayer graphene. The arrow indicates the direction of an applied external electric field. 

therefore essential. 

Here we systematieally investigate the effeet of physisorbed hexaeene and perfluorohexaeene 
(fluorinated hexaeene) on the eleetronie properties of bilayer graphene using first prineiples den¬ 
sity functional theory simulations. We use perfluorohexacene as an effective tool to down-shift the 
molecular energy levels relative to hexacene, and induce significant n-n* interactions and symme¬ 
try breaking in bilayer graphene. We examine how the functional groups and adsorption geometry 
of molecules influence the stability and the electronic structure of the bilayer graphene-molecule 
system. We show that the adsorption of hexacene and perfluorohexacene on bilayer graphene leads 
to a significant charge redistribution and the formation of localized states in graphene. By applying 
an external electric field to bilayer graphene adsorbed with hexacene and perfluorohexacene, we 
demonstrate that the induced localized states in graphene can be effectively controlled, potentially 
providing a new strategy for graphene-based sensors for a selective sensing of weakly adsorbed 
molecules. 
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Methods 


To obtain optimized geometries and the eleetronie struetures of all graphene-moleeule systems 
eonsidered in our study, we employed first prineiples density funetional theory as implemented in 
Vienna Ab Initio Simulation Paekage. We used the generalized gradient approximation of the 
Perdew-Burke-Emzerhof form for the eleetron exehange-eorrelation funetional.^* The eore and 
valenee eleetrons were treated using projeetor augmented wave (PAW) seheme*^ with a kinetic 
energy cut-off of 600 eV for the plane-wave basis set. Since the generalized gradient approxima¬ 
tion does not fully account for long-range dispersion interactions,"^® we used a Grimme’s semi- 
empirical functional to account for these interactions in the weakly bound graphene-molecule 
system. To benchmark the accuracy of this functional, we obtained the the equilibrium interlayer 
distance for pristine bilayer graphene to be 3.23 A, which is in good agreement with the experi¬ 
mental values. We used a periodic 8x4 graphene supercell for investigating the adsorption of 
hexacene and perfluorohexacene on Bernal-stacked bilayer graphene. A single molecule in this su¬ 
percell represents a nearly monolayer coverage for hexacene and perfluorohexacene on graphene 
with a molecular density of 9.846 x 10'^^ mol/cm^. This magnitude of molecular density is repre¬ 
sentative of the reported coverage of the aromatic molecules deposited on graphitic or graphene 
surfaces in various experimental studies. In each case, the periodic images were separated by 
a 30 A vacuum, which was found to be large enough to avoid any spurious interactions between 
the periodic images. All structures were fully relaxed until the ionic forces were smaller than 0.01 
eV/A. Gaussian smearing was used for geometry relaxations, while Bloch tetrahedral smearing 
was employed for subsequent calculations of electronic structures.^® Finally, for accurate calcula¬ 
tions of the electronic structures, we used a fine 6x12x1 F-centred grid for sampling the Brillouin 
zone. 

First, we individually examined the adsorption of hexacene and perfluorohexacene on bilayer 
graphene using a single-molecular adsorption configuration as shown in Fig. 1(a). We considered 
two stacking sequences for this configuration, namely, ABA and ABC, using the same notation as 
in the case of few-layer graphene. We also considered the simultaneous adsorption of hexacene 
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and perfluorohexacene in a dual-molecular configuration shown in Fig. 1(b). Finally, the influence 
of external electric fields normal to the plane of graphene was investigated by introducing dipolar 
sheets at the center of supercells. 

Results and Discussion 

Electronic Structure of Hexacene and Perfluorohexacene 

We first obtained the optimized geometries of hexacene and perfluorohexacene molecules and 
calculated their electronic structure. The relaxed geometries of hexacene and perfluorohexacene 
are shown in Fig. 2(a). The average C-H bond length for hexacene and the average C-F bond length 
for perfluorohexacene were calculated to be 1.09 A and 1.35 A respectively. These bond lengths 
are in close agreement with previous structural calculations reported by Kadantsev et al. Figure 
2(b) depicts an energy level diagram showing the highest occupied molecular orbital (HOMO) 
and the lowest unoccupied molecular orbital (LUMO) of hexacene and perfluorohexacene with 
respect to the position of the Fermi level of graphene. The HOMO-LUMO band gaps of hexacene 
and perfluorohexacene are narrow, namely, 0.80 eV and 0.59 eV, respectively. The calculated 
electronic band gap of these molecules is smaller than the experimentally determined band gaps 
by 0.5 eV - 0.6 eV due to a systematic underestimation of band gap values of semiconducting 
materials obtained by DFT calculations employing the generalized gradient approximation. 36.49,50 

Single- and Dual-Molecular Adsorption on Bilayer Graphene 

Next we obtained optimized ground state configurations of hexacene and pefluorohexacene ad¬ 
sorbed on bilayer graphene in single as well as dual-molecular configurations shown in Fig. 1. 
In each case, the adsorption energy is calculated as AE = Egraphene/moiecuie -Egraphene - ^molecule, 
where Egraphene/moiecuie is the total energy of the fully-relaxed graphene-molecule supercell, while 
Egraphene ^nd ^molecule is the energy of bilayer graphene and isolated molecules in the same su¬ 
percell, respectively. A negative value of the adsorption energy indicates an exothermic, thermo¬ 
dynamically favorable adsorption. The values for adsorption energy and the average adsorption 
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Figure 2: (a) Relaxed geometries of hexacene and perfluorohexacene. The C-C, C-H and C-F bond lengths 
(in A) are indicated, (b) The energy level diagram of the band alignment between the Fermi level of bilayer 
graphene, and the HOMO and LUMO of hexacene and perfluorohexacene (in eV). The energy of vacuum is 
regarded as zero. 

distance from the nearest graphene layer for eaeh eonfiguration are summarized in Table 1 . We 
find that the adsorption of both hexaeene and pefluorohexaeene on bilayer graphene is thermo- 
dynamieally favorable in all eonfigurations considered. The resulting adsorption distanee of the 
moleeules from graphene is elose to the interlayer distanee between graphene layers (~3.23 A), 
indieating that graphene and the moleeules are bound by weak van der Waals dispersion interae- 
tions. The adsorption distanee and the adsorption energies are also nearly unaffeeted by the ehange 
of stacking sequence from ABC to ABA, also suggesting that the moleeule-bilayer graphene inter¬ 
actions are largely confined between the molecule and the adjaeent graphene layer. 

Our results show that pefluorohexaeene binds more strongly to bilayer graphene than hexae¬ 
ene (adsorption energy of -2.25 eV vs. -1.76 eV) in a single-moleeular adsorption eonfiguration. 
The stronger binding of pefluorohexaeene eorrelates with a higher eleetron affinity and ehemieal 
reactivity of pefluorohexaeene relative to hexaeene. Finally, in the ease of simultaneous, dual- 
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Table 1: Adsorption energy and the average adsorption distance for monolayer coverage of hexacene 
and perfluorohexacene on bilayer graphene in single- and dual-molecular adsorption configurations. 



Stacking 

Adsorption 

Average adsorption distance (A) 


sequence 

energy (eV) 


Hexacene 

ABA 

-1.756 

3.19 

ABC 

-1.751 

3.19 

Perfluorohexacene 

ABA 

ABC 

-2.246 

-2.239 

3.16 

3.16 

Dual-molecular adosprion 

ABAB 

-4.015 

3.20 (hexacene) 

3.17 (perfluorohexacene) 

(hexacene-i-perfluoroxacene) 

ABCA 

-4.005 

3.19 (hexacene) 

3.16 (perfluorohexacene) 


molecular adsorption of hexacene and perfluorohexacene on bilayer graphene, the adsorption en¬ 
ergy of the total hexacene and perfluorohexacene complex is nearly equal to the sum of the adsorp¬ 
tion energies of individual molecules in single-molecular adsorption. This observation suggests 
that the adsorption system can be relatively easily adjusted from a single-molecular configuration 
to a dual-molecular configuration, or vice versa, at a very low energy cost. The adsorption dis¬ 
tances of hexacene and perfluorohexacene in dual- molecular adsorption are found to be slightly 
larger than for the single-molecular adsorption cases. 


Electronic Structure of Bilayer Graphene upon Adsorption 

Following the structural optimisation, we next investigate the electronic properties of bilayer graphene 
upon molecular adsorption. Our results show that the molecular stacking sequence did not have 
any effect on the density of states curves. Consequently, Fig. 3 presents the partial density of states 
(PDOS) for bilayer graphene with single-molecular adsorption with ABA stacking and for dual- 
molecular adsorption with ABCA stacking sequence. We find that the physisorbed hexacene and 
perflurorohexacene has a negligible effect on the n states of bilayer graphene in the vicinity of 
the Dirac point. No shift in the Fermi level of bilayer graphene was observed upon adsorption of 
hexacene or perfluorohexacene. The magnitude of the band gap induced in bilayer graphene upon 
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Figure 3: (a) Total density of states of a pristine bilayer graphene, (b-d) Partial densities of states (PDOS) 
of bilayer graphene upon single-molecular adsorption of hexacene (b), single-molecular adsorption of per- 
fluorohexacene (c), and dual-molecular adsorption of hexacene and perfluorohexacene (d). In all cases, the 
molecular concentration is 9.85x10'^^ mol/cm^ and represents monolayer coverage. (Hex)H, (Hex)L, and 
(P-Hex)H, (P-Hex)L denote the localized states induced by hybridization with HOMO and LUMO of hexac¬ 
ene, and HOMO and LUMO of perfluorohexacene, respectively. The Fermi level is set to zero in each case, 
(e-f) Partial charge density plots for the localized states in PDOS of bilayer graphene upon single-molecular 
adsorption of hexacene (e) and single-molecular adsorption of perfluorohexacene (f). The isosurface level 
is set to 0.0003 e/A^. 


adsorption of these moleeules is also negligible (3 meV and 5 meV for hexaeene and perfluroro- 
hexaeene, respectively). These observations point to a weak interaction between graphene and the 
molecules. Nevertheless, the low lying highest occupied and lowest unoccupied molecular orbitals 
(HOMO and LUMO) of adsorbed molecules hybridize with %!%* states of graphene, giving rise 
to two localized states near the Fermi level of bilayer graphene (see Fig. 3(a-d)). These states 
were found in PDOS of bilayer graphene at -0.18 eV and 0.61 eV after the adsorption of hexacene 
and at -0.43 eV and 0.18 eV after the adsorption of perfluorohexacene. The position of HOMO 
and LUMO of hexacene and perfluorohexacene states is slightly altered in comparison to that of 
isolated molecules (refer Fig. 2 (b)). To gain further insight into the localized states in graphene. 
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we also plotted the corresponding partial charge densities by integrating the charge density in an 
energy range ±0.02 eV around the localized peaks. The partial charge density plots (Fig. 3(e,f)) 
show the shape of HOMO/LUMO of hexacene/perfluohexacene, and present the signature of hy¬ 
bridization between graphene and the adsorbed molecules. It is evident that the induced localized 
states in graphene are located on the nearest carbon atoms of the top graphene layer. The pres¬ 
ence of these hybridized states near the Dirac point of graphene implies that both electrons and 
holes can be injected from graphene to molecules at a relative low energy cost. In dual-molecular 
adsorption configuration (Fig. 3 (d)), the PDOS can be regarded as a superposition of the energy 
states from the single-molecular adsorption on graphene. Moreover, we observe that a 8 meV band 
gap is opened, equal to the sum of the band gap values of individual molecules. This indicates 
that the interaction between bilayer graphene and acene molecules in dual-molecular configuration 
is essentially governed by the interaction of a single graphene layer with the adjacent adsorbed 
molecule. 

In order to further assess the influence of molecular adsorption on the electronic properties 
of graphene, we calculated the charge density difference as defined by Ap = Pgraphene/moiecuie — 
Pgraphene Pmoleculej where pgraphene/moleculej Pgraphene Pmolecule ^6 the electronic charge densi¬ 
ties of the adsorbed system, isolated graphene and the molecule, respectively. With this definition, 
a positive value of Ap indicates an accumulation of electronic charge and a negative value indi¬ 
cates a charge depletion. The distribution of the charge density difference for single-molecular 
adsorption of hexacene and perfluorohexacene is shown in Fig. 4. In the case of adsorption of 
hexacene on bilayer graphene (Fig. 4 (a)), the charges are depleted from the region 0.6 A-1.0 
A above the top graphene layer and accumulated close to hexacene in the region 2.6 A- 2.8 A 
above the graphene layer. This charge redistribution primarily arises from the electrostatic inter¬ 
action between the aromatic rings of hexacene and graphene—the interaction between hydrogen 
atoms of hexacene and carbon atoms of graphene or hexacene is found to be negligible. Overall, 
the interaction between hexacene and bilayer graphene is not strong enough to lead to a signifi¬ 
cant direct charge transfer between graphene and hexacene, but only to a charge redistribution on 
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Figure 4: The distribution of charge density difference for single-molecular adsorption of (a) hexacene and 
(b) perfluorohexacene on bilayer graphene at monolayer coverage. Red and green isorufaces indicate the 
accumulation and depletion of electrons at a level of 0.0003 e/A^, respectively. 


carbon atoms of the molecule and the nearest graphene layer. In contrast to hexacene, the charge 


redistribution is significantly different for perflurorohexacene physisorbed on bilayer graphene as 
shown in Fig. 4 (b). Since fluorine atoms are strong electron-attracting groups, the n electrons are 
polarized away from the aromatic rings leading to a relatively electron-deficient aromatic core of 
perfluorohexacene. Similarly, a strong interaction between fluorine and graphene gives rise to a 
significant charge depletion from the carbon atoms in graphene close to fluorine atoms in perflu¬ 
orohexacene, indicating a net charge transfer from graphene to perfluorohexacene. By comparing 
the charge distribution of hexacene and perfluorohexacene on bilayer graphene, it is evident that in¬ 
teraction between graphene and perfluorohexacene is largely controlled by the presence of fluorine 
functional groups. 

To quantify the charge transferred from graphene to adosprbed molecules, we plotted the varia¬ 
tion of the planar-averaged and the integrated planar-averaged charge density difference as a func¬ 
tion of the distance from the basal plane of graphene as shown in Fig. 5. The planar-averaged 
charge density difference Apavg(z) along a plane parallel to the basal plane of graphene is obtained 
by integrating the charge density difference across the plane, whereas the integrated charge den¬ 
sity difference Ap(z) is simply calculated by integrating Apavg(z) curve from the boundary of the 
periodic box to the position of the plane. The extremum in the integrated charge density difference 
curve in the region between graphene and the adsorbed molecule (indicated by red lines in Fig. 5) 
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Figure 5: The planar-averaged (Apavg(z)) and integrated planar-averaged (Ap(z)) charge density difference 
curves for single-molecular adsorption of (a) hexacene and (b) perfluorohexacene on bilayer graphene at 
monolayer coverage. The positions of the top graphene layer (TG) and the molecule are marked by green 
and yellow dashed lines, respectively. Blue solid lines denote the neutral plane. 

denotes the neutral plane and the magnitude of the net eharge transfer. Aeeording to this analysis, 
graphene donates 0.035 per moleeule to hexaeene and a larger fraction 0.050 per molecule 
to perfiurohexacene. The magnitude of net charge transfer obtained here is qualitatively consistent 
with Bader charge population analysis, which yields a net charge transfer of 0.02 and 0.18 
e~ per molecule from graphene to hexacene and perfluorohexacene, respectively. The differences 
in the magnitudes of the net charge transfer obtained by these two methods can be attributed to 
the different ways of calculating the net charges on each ion. Bader charge analysis takes into 
account both core charges and valence charges, while the charge density difference calculation is 
more suitable for interpreting the charge redistribution close to the Fermi level. Nevertheless, the 
small magnitude of charge transfer for both molecules is indicative of a weak interaction between 
the molecules and bilayer graphene. 
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Finally, few recent experimental studies have reported that as the density of the adsorbed acene 
molecules is increased beyond the near-monolayer coverage considered here, the molecules tend 
to tilt rather than maintain a planar orientation on graphene/^’^^ In order to confirm this observa¬ 
tion, we have also studied the adorption of hexacene on bilayer graphene at a nominal molecular 
concentration of 1.31x10'^® mol/cm^ using a 8 x 3 graphene supercell. Figure SI in the Sup¬ 
porting Information shows the adsorption geometry as well as the electronic interaction between 
hexacene and bilayer graphene at this coverage. We find that hexacene shows a remarkable 11° 
tilt with respect to the basal plane of graphene, in qualitative agreement with the experimental 
observations.This tilt arises due to a stronger repulsive intermolecular interaction between 
neighbouring hexacene molecules at high coverage. The cofacial n-n interactions that contribute 
to the stability of hexacene adsorption at lower concentrations are disrupted, causing the adsorp¬ 
tion at high coverage to be less energetically favorable (adsorption energy -1.534 eV per molecule). 
Moreover, the electronic interactions between graphene and hexacene also vary spatially, leading 
to an asymmetric charge redistribution pattern as shown in Fig. S1. Compared to the correspond¬ 
ing patterns at low coverages, a significantly larger charge rearrangement is observed in the region 
of the molecule where hexacene is closer to the graphene than in the region where hexacene is 
away from the graphene. This imbalance in the charge redistribution breaks the local symmetry of 
bilayer graphene, inducing a 54 meV band gap. 

Effect of the Applied External Electric Field 

The results presented in earlier sections show that in general, the adsorbed aromatic acene molecules 
interact weakly with bilayer graphene, leading to the formation of localized states and a weak 
p-type doping of graphene. Since external electric fields can enhance the interactions between 
adsorbate and substrate, next we investigate whether the external electric field could be ef¬ 

fectively utilized to tune the electronic structure and molecule-specific localized states in bilayer 
graphene. Specifically, we studied the effect of the electric field in the range of -3 to 3 eV/nm, ap¬ 
plied perpendicular to the basal plane of bilayer graphene in the dual-molecular hexacene/bilayer 
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Figure 6: Effect of external electric field on the electronic structure of bilayer graphene in dual-molecular 
adsorption of hexacene and perfluorohexacene at monolayer coverage, (a) PDOS of bilayer graphene as 
a function of external electrical field (in eV/nm). The peaks induced by hybridization with HOMO and 
LUMO states of hexacene and perfluorohexacene are denoted. The grey region indicates the band gap 
induced in bilayer graphene. The Fermi level is set to zero, (b) Positions of the localized states in bilayer 
graphene induced by the hybridization with HOMO and LUMO states of hexacene and perfluorohexacene. 
(c) Induced charges in hexacene (Hex), perfluorohexacene (P-Hex), and top and bottom graphene layers 
(TG, BG), calculated by Bader charge population analysis. 


graphene/perfluorohexacene adsorption configuration at a monolayer coverage as shown in Fig. 
1(b). In our notation, a positive electric field is oriented towards hexacene from perfluorohexac¬ 
ene. Figure 6(a) presents PDOS of bilayer graphene as a function of the strength of the field. The 
application of the electric field leads to the opening of a considerable band gap in bilayer graphene, 
as well as to the shift in the energy levels of the localized states arising due to the hybridization 
with HOMO and LUMO states of the molecules. The external electric field causes an accumulation 
of electrons in one layer and a depletion of electrons in the other layer of bilayer graphene, thus 
breaking the interlayer symmetry. The band gap generated in bilayer graphene is thus a result of 
the interplay between the field-induced interlayer symmetry breaking and the asymmetric charge 
transfer between graphene and the adsorbed hexacene and perfluorohexacene molecules. In several 
cases, the magnitude of band gap is difficult to determine from the PDOS of bilayer graphene alone 
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as the gap region is oeeupied by the loealized states. In order to eorreetly identify the band gaps, 
we have ealeulated the eleetronie band strueture for eaeh ease as shown in Fig. 7. The loealized 
states ean then be readily identified from these band strueture diagrams as the flat bands between 
the n and n* bands of bilayer graphene. 

It is evident that both the shape of the PDOS as well as the magnitude of the band gap in bilayer 
PDOS shows a strong dependenee on the strength and the direetion of the external eleetrie field. 
For example, the band gap inereases linearly with the eleetrie field till a maximum of ~250 meV 
when the field is oriented from hexaeene towards perfluorohexaeene (that is, a negative eleetrie 
field). For a positive field, the band gap reaehes ~200 meV for the field magnitude of 1 eV/nm, 
then reduees to ~115 meV for the field of 2 eV/nm or greater. This observed eontrast in the trend 
between negative and positive fields ean be attributed to the distinet eharge transfer behavior of 
the adsorbed hexaeene and perfluorohexaeene moleeules as shown in Fig. 6 (b). For instanee, for 
fields E > 2 eV/nm, the HOMO of perfluorohexaeene and the LUMO of hexaeene are pinned in 
the vieinity of the Fermi level. The pinning of these loealized states near the Fermi level enhanees 
the driving foree for transferring eharge from hexaeene to perfluorohexaeene, whieh is favored by 
the applieation of a positive eleetrie field. Therefore the eharge transfer between hexaeene and 
perfluorohexaeene is maximized as shown in Fig. 6 (e). The eharge inequivalenee between two 
graphene layers is redueed for positive eleetrie fields greater than 1 eV/nm, leading to a saturation 
of the band gap. For negative fields on the other hand, the eharge transfer trend is reversed. The 
energy differenee between the HOMO of perfluorohexaeene and the LUMO of hexaeene inereases 
with the strength of the negative eleetrie field. Therefore the eharge transfer between the two 
graphene layers is less affeeted and the magnitude of band gap inereases rapidly with the magnitude 
of the field. 

Figure 7(a) also shows a eomparison between the eleetronie band strueture diagrams of bi¬ 
layer graphene adsorbed with hexaeene and perfluorohexaeene to those with the pristine bilayer 
graphene. Our obtained band struetures for pristine graphene are in good agreement with previous 
works.It is evident that due to the deformation of n and Tt* bands, the field-indueed band gap 
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Figure 7: (a) Electronic band structure diagrams for dual-molecular adsorption of hexacene and perfluo- 
rohexacene with monolayer coverage on bilayer graphene as a function of the strength of the electric field 
(in eV/nm). The corresponding diagrams for pristine bilayer graphene are also shown for comparison. The 
TT and 71* bands of bilayer graphene are labeled in red color. The inset shows the first Brilloun zone of the 
supercell compared against that of the unit cell of graphene, (b) The variation of the band gap at K (Ek, 
blue) and the true band gap along EX line in the Brillouin zone (Ex, yellow) with the strength of the electric 
field. 
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in pristine bilayer graphene is no longer loeated at the K point, but instead along F-X line of the 
Brillouin zone. However, the n and n* bands of bilayer graphene with dual-moleeular adsorption 
are less deformed elose to the Fermi level due to the sereening by moleeules, indieating that the 
energy dispersion relationship in the vieinity of the K point is relatively well preserved. This ean 
be elearly seen by eomparing the magnitude of band gap at K (Ek) and the true band gap along 
F-X line (Ex) presented in Eig. 7(b). It ean be seen that in general, the range of band gaps that ean 
be indueed in bilayer graphene with moleeular adsorption (100 meV - 250 meV) is eomparable to 
the pristine bilayer graphene (200 meV - 300 meV). Overall, the band struetures follow a similar 
trend as the DOS plots of dual moleeular adsorption on bilayer graphene shown in Eig. 6 (a), where 
graphene n and n* states hybridize with the loealized HOMO and EUMO bands of the moleeules, 
forming new loealized bands with a flat dispersion. These results show that the energy of these 
loealized states ean be varied as a funetion of the strength and polarity of the applied eleetrie 
field, thereby modifying the eleetronie strueture of bilayer graphene. Einally, for eompleteness, 
we briefly eompare the eleetronie strueture of bilayer graphene with dual-moleeular adsorption 
under the eleetrie field to that of monolayer graphene. Eigure S2 in the Supporting Information 
shows band strueture diagrams for dual-moleeular adsorption of hexaeene and perfluorohexaeene 
on monolayer graphene as a funetion of the strength of the eleetrie field. We find that the eleetronie 
strueture of monolayer graphene with dual-moleeular adsorption is signifieantly different from the 
bilayer graphene for negative eleetrie fields. The loealized states indueed in monolayer graphene 
due to hybidization with the HOMO/EUMO states of the adsorbed moleeules are relatively further 
away from the Eermi level of graphene. Moreover, negative eleetrie fields lead to mueh lower band 
gaps in monolayer graphene (less than 50 meV). Eor positive eleetrie fields, both monolayer and 
bilayer graphene show a similar magnitude of the band gaps (~120 meV). When eompared with 
Eig. 7, these observations highlight the interplay between the eharge transfer and the breaking 
of the symmetry between the top and bottom layer in bilayer graphene. In general, these results 
suggest that by the applieation of external eleetrie field to bilayer graphene with dual moleeular 
gating, the eleetronie strueture ean be more flexibly eontrolled, leading to opening of eonsiderable 
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band gaps. 


Conclusions 

Here we have used self-consistent density functional theory calculations to study the effect of 
physisorption of hexacene and its fluorine derivative, perfluorohexacene, on the electronic struc¬ 
ture of bilayer graphene. We find that although the overall interaction between graphene and 
molecules is weak, the adsorption of these molecules results in a significant charge redistribution. 
This charge redistribution gives rise to the hybridization of HOMO/LUMO energy levels of the 
molecules with the n electrons of graphene, leading to the formation of localized states in bilayer 
graphene. We have shown that the external electric fields can be used to tune the electronic prop¬ 
erties of graphene-molecule system, effectively opening large band gaps of the order of 250 meV 
in bilayer graphene. Furthermore, external electric fields can also infleunce the energies of the 
localized states in graphene, an effect that can be utilized in organic field effect transistor (OFET) 
devices by aligning the electronic states of acene channels with that of graphene electrodes. This 
effect can also be potentially useful in the sensing of different organic molecules on the surface 
of graphene transistors. Graphene transistors have proven to be extremely sensitive sensors,but 
their selectivity remains a major problem for their practical use. In summary, we have shown that 
hexacene, a stable and high mobility organic electronic material, and its derivatives, are promising 
candidates for surface electronic structure modification of graphene for potential applications in 
organic electronics and sensing. 
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Supporting Information 

Additional figure for the adsorption of hexacene on bilayer graphene in the 8x3 graphene supercell 
and band structure diagrams for dual-molecular adsorption of hexacene and perfluorohexacene on 
monolayer graphene. 
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c 




Figure SI. Single-molecular adsorption of hexacene on bilayer graphene at a high 
(LSlxlO'*** mol/cm^) molecular concentration, (a) Adsorption geometry with 11“ tilt 
angle, (b) PDOS of bilayer graphene. The grey region indicates a 54 meV band gap, 
the Fermi level is set to zero, (c) Charge density difference along a plane parallel to 
the basal plane of graphene. The position of the cutplane is 0.62 A above the top 
graphene layer. The color bar indicates the level of isosurface (unit in 10'^ e/A^). 
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Figure S2. Electronic band structure diagrams for dual-molecular adsorption of 
hexacene and perfluoroliexacene with monolayer coverage on monolayer graphene as 
a function of the strength of the electric field (in eV/nm). The n and n* bands of 
bilayer graphene are labelled in red color. 










































